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ABSTRACT 

We extend our previous study of the stellar population of L1641, the lower-density star-forming 
region of the Orion A cloud south of the dense Orion Nebula Cluster (ONC), with the goal of testing 
whether there is a statistically significant deficiency of high-mass stars in low-density regions. Previ- 
ously, we compared the observed ratio of low-mass stars to high-mass stars with theoretical models 
of the stellar initial mass function (IMF) to infer a deficiency of the highest- mass stars in L1641. We 
expand our population study to identify the intermediate mass (late B to G) L1641 members in an 
attempt to make a more direct comparison with the mass function of the nearby ONC. The spec- 
tral type distribut i on and the K-band luminosity function of L1641 are similar to those of the ONC 
(jHillenbrandl 119971: iMuench et al.l I20Q2D , but problems of incompleteness and contamination prevent 
us from making a detailed test for differences. We limit our analysis to statistical tests of the ratio of 
high-mass to low-mass stars, which indicate a probability of only 3% that the ONC and the southern 
region of L1641 were drawn from the same population, supporting the hypothesis that the upper mass 
end of the IMF is dependent on environmental density. 

Subject headings: surveys — stars: formation — stars: luminosity function, mass function — stars: 
pre-main sequence 



1. INTRODUCTION 

Does the formation of high-mass stars depend system- 
atically on the environmental density of the star-forming 
cloud? In a r ecent review of this important question, 
iBastian et al.l (|201Cll ) found no clear evidence for varia- 
tions of the stellar initial mass function (IMF) as a func- 
tion of initial conditions, but also concluded that further 
study in specific local and extragalactic environments is 
warranted. Even if the IMF is reasonably stable on a 
galactic scale, as a result of averaging over a variety of 
environments, it would still be of interest to discern IMF 
variations in specific regions as a clue to the processes of 
star formation. 

Studies of very young regions are necessary to make a 
direct connection between environmental conditions and 
the production of massive stars. The most detailed in- 
vestigations of the most populated ne arby you ng clusters 
- the Orion Nebula Cluster fONO (IHillenbr and 1997; 
Muench e t al][200l iDaRio et al.ll2010D and NGC 2264 
Sung et all I2004l " are consistent with the the canonical 
IMF in general and with the Sal peter slope i n the upper- 
mass end in particular (Bastian et al.|[20Tol ). Producing 
a statistically-significant test of the high-mass IMF is 
much more difficult in low- density regions (~ 1 - 10 
stars pc^^), such as Taurus or Chamaeleon, simply be- 
cause there are too few low-mass stars to make a strong 
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test (e.g. ILuhman et al.l [20091 for Taurus and ILuhmanI 
120071 for Chamaeleon I). Therefore, we need a survey of 
a large, well-populated low-density region with spectro- 
scopic confirmation of membership to rule out line-of- 
sight contamination. It is also important to compare the 
high- and low-density populations directly rather than 
comparing with inferred "universal" IMF forms, as the 
results may depend on the IMF form used and the statis- 
tical significance of d isagreement can be difficult to infer 
(Bas tian et al.ll2010t) . 

The Orion A molecular cloud is an ideal site to study 
star-formation in both high and low-density environ- 
ments. While the ONC is one of the most thoroughly 
observed young dense cluster (the stellar surface den- 
sity is ~ 1000 pc^^ for the Trapezium cluster and ~ 200 
pc~^ for the whole ONC), the outlying regions provide 
us with a substantial low-density population. L1641, 
the low-density star- forming region south of the ONC, 
has a stella r population com parable to the ONC in size 
{N > 1000: iHsu et~alll2012t hereafter, Paper I), but has 
no large scale clustering and a much lower stellar surface 
density of ^ 10 pc~^, which makes it the best low-density 
region to study the high-mass IMF. It also has the ad- 
vantage of being at approximately the same' distance a s 
the ONC (~ 414Dc: lMmren et al.ll2007l : IKim et al.l[2()?)8h . 
allowing the two populations to be compared directly. 

In Paper I we presented an optical spectroscopic and 
photometric survey of the low-mass population in L1641 
to test whether the previously-recognized lack of high- 
mass stars is statistically significant. Co mbining the op- 
tical d ata to the Spitzer/IRAC survey bv lMegeath et al.l 
(|2012D , we identified and spectral- typed nearly 900 mem- 
bers; as we are unable to observe the heavily-extincted 
members, we estimated that L1641 may contain as many 
as ^1600 stars down to O.IM0. Based on the large 



2 



Hsu et al. 



number of low-mass stars, we concluded that the lack 
of stars e a rlier t han B4 is inconsis tent with the canonical 
IChabried ([200l or lKroupal (l200l IMFs to 3-4 a signifi- 
cance. 

In this paper we expand our optical photometric and 
spectroscopic sample to the intermediate mass (late B to 
G) members to improve our estimates of the mass func- 
tion. We first describe our observations and data reduc- 
tion in § [2l Then, in § |3l we describe the observational 
results and attempts to use proper motions and radial ve- 
locities to reduce foreground/background contamination. 
In addition, we estimate the age of L1641 to be about 3 
Myrs (see l3.6p . so late O and early B stars should still be 
present. In § |4l we compare the spectral type distribu- 
tion and K-band luminosity function (KLF) of L1641 to 
that of the ONC, finding that we cannot make detailed 
tests for differences between the two regions given prob- 
lems of incompleteness and contamination. In § [Sj we 
relate our results to t he mass of most massive star and 
cluster mass relation (jWeidner fc Kroupall2006D and dis- 
cuss the challenges in searching for density dependence 
of the IMF. Our results are then summarized in § [6l 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Target Selection 

We aimed to target all the high to intermediate mass 
stars in L1641 with moderate extinction, regardless of 
whether they have disks. We selected o ur targets based 
on their 2MASS colors and USNOb ( Monet et all [200l 
photometry because our photometric survey (described 
in Paperl) saturates at around 1^-^12, which corresponds 
to late G to early K stars in LI 641 in less extincted re- 
gions. The targets were selected using the following pro- 
cess: 

1. We plotted the J vs. J-H color- magnitude diagram 
(CMD) of ah objects with 2MASS magnitudes in 
L1641 (top panel of Figure[T]). We then selected the 
targets that lie above the extinction vector orig- 
inating from the location of a IM© star on the 
iBaraffc et al. (1998) 3 Myrs isochrone. This sam- 
ple included objects, if at Orion's distance, that are 
more massive than IM©. 

2. We excluded stars that already had optical spectral 
types from Paper I, which left us with 446 stars. 
In the top panel of Figure [U the blue circles show 
objects with spectral types from Paper I. 

3. We then compared the coordinates of our targets to 
the USNOb catalog. If there is not a match, then 
the star is too faint for our optical spectroscopy. 
We selected only objects that have B-R < 2.5 and 
B brighter than 15 in the USNOb catalog for spec- 
troscopic foUowup. 

A total of 136 potential intermediate-mass members 
were observed spectroscopically p.2p . Their positions 
on the J vs. J-H CMD are shown by the stars in the 
top panel of Figure [TJ The bottom panel of Figure [T] 
shows the combined sample of stars observed in this pa- 
per with the Ohio State Multi-Ob ject Spectrograp h (OS- 
MOS) on the MDM 2.4m (jStoU et al...2010; .Martini et al.l 
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Figure 1. Top: J vs. J-H color-magnitude diagram (CMD) used 
to select targets from the 2MASS catalog. The grey shaded region 
indicates our initial cut for selecting targets. The blue solid circles 
are objects with spectral types from Paper I. The purple stars 
are objects that we obtained spectra with OSMOS in this study. 
Bottom: USNOb B vs. B-R CMD of targets selected in the top 
panel. The solid triangles are targets observed in this work, the 
open triangles are targets observed in Paper I and the the black 
dots are objects that were not observed. 

120111 ) and stars from Paper I. The solid triangles are tar- 
gets observed in this work, the open triangles are tar- 
gets observed in Paper I. There are also a small num- 
ber of stars that were not observed (black dots). These 
targets are too faint for the OSMOS spectroscopic ob- 
servations. They were not observed with Hectospec or 
IMACS in Paper I because they either did not satisfy 
the selection criteria based on their position on the V vs. 
V-I CMD, including stars that we do not have I band 
photometry due to saturation. Also shown in this fig- 
ure are the exti nction vector corresp o nding to Av=l and 
the ZAMS from'Leieune & Schaerer C 2001D (with models 
from [Schaller et al.,,1992, ). Figure [2] shows the positions 
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Figure 2. Positions of all stars observed witii OSMOS spec- 
troscopy, overlaid on ^^CO map from Bally et al. (1987). The blue 
boundary lines repres ent the fie lds covered by the Spitzer/IKAC 
survey of Orion ( Meg eath et al.| [2012). The circles showr the fields 
of the OSMOS photometry from Paper I. 

of the observed stars on the sky. 

2.2. Optical Spectroscopy of Intermediate Mass Stars 

We obtained optical spectroscopy of selected objects 
(i )2.ip with OSMOS. The observations were taken with 
center slit on the VPH grism, which yields a spectral 
range of 3100A to 5950A and efficiency peak at 5000A. 
The 1.4" slit was used, giving a spectral resolution of 
R=1400, or r^2 - 4A for the spectral range. 

The data is first processed with the idl program proc4k 
written by Jason Eastman to subtract overscan and then 
bias subtracted, flat-field corrected with the MIS fiat 
lamp to remove small-scale variation of pixel responses. 
Xenon comparison lamp (exposure time of 250s) is used 
for wavelength calibration as recommended for the cen- 
ter slit. Then the IRAF apall program is used to trace 
aperture and extract the spectra. 

2.3. Optical Photometry of Intermediate Mass Stars 

We used OSMOS in photometric mode to obtain 
BVRIc band photometry in the Johnson-Cousin system. 
The photometry observations were all done on the night 
of Jan 1 - Jan 2, 2012. The conditions were photometric 
with seeing ^ 1.2" in the middle of the night. 

Since the stars we are interested in are very bright stars 
that saturate in our previous photometric observations, 
we need to use an even shorter exposure time. (The 
saturation limit is 1^12 for 5s exposures; see Paper I for 
details.) We therefore used exposure time of 1 second for 
all four bands and took another 5 second exposure for the 
B band image. We observed the targets once with the 
telescope in focus and then once with the telescope out of 
focus and the images of stars are shaped like doughnuts. 



This strategy allows us to obtain photometry of brighter 
stars without saturating. The FWHM of the in focus 
images are around 1.1 - 1.4", depending on the seeing, 
and the out of focus images have FWHM of 2 - 2.5 ". The 
radius for the aperture photometry is 9 pixels, or 5". To 
minimize errors due to shutter timing, we place the stars 
near the center of the field and use only the center Ik 
X Ik of the 4k x 4k CCD. The difference between the 
center and the edge of the field is about 2% in this Ik x 
Ik field. The 2x2 binning was used, which gives a plate 
scale of 0.55" per pixel. For calibrations, we observed 
Landoh fields SA92, SA95, SA98 and SAlOl throughout 
the night with the fuU field of view of OSMOS. 

Each CCD frame was first corrected by overscan us- 
ing the IDL program proc4k written by Jason Eastman. 
Basic reduction was performed following the standard 
procedure using IRAF. Sky fiats were used in the flat- 
field correction. We then obtained aperture photometry 
with the IRAF phot package and manually choose our 
target stars or the Landolt standard stars. If the 5 sec- 
ond B band exposure is saturated, we use the magnitude 
obtained from the 1 second exposure. The photometric 
errors calculated by IRAF phot package are smaller than 
0.04 mag in all cases. The rms departures of the Landolt 
stars from the calibration equations are ^ 0.03 mag for 
all bands. We therefore expect our photometric errors to 
be less than 0.05 mag. 

We were unable to obtain optical photometry for stars 
brighter than 9th magnitude in V due to saturation. We 
therefor e use the B fc V photometry from the Tycho-2 
catalog (|H0g et all [20000 . The Tycho-2 catalog is 99% 
complete to V 11.0 and the magnitude error is 0.013mag 
for stars brighter than 9th magnitude and O.lmag for 
all other stars. The Tycho-2 photometry is transformed 
to the Johnson system by using the following equation 
(jPerrvman & ESAl[T99l : 

Vj = Vt- 0.090(E - V)t 

{B ~V)j = 0.850{B -V)t (1) 

2.4. Echelle Observations of Selected F & G Stars 

To better distinguish foreground F & G stars from 
members (see Sec 13. 3p . we observed 20 of the F & G 
stars with MIKE, the Magellan Inamori Kyocera Echelle 
on the Magellan Clay telescope on Feb 4 and 5, 2012. 
We used both the blue side (wavelength coverage: 3200 

- 5000A) and the red side (wavelength coverage: 4900 

- lOOOOA) simultaneously. The observations were taken 
with 0.7" slits, which gives a resolving power of 42000 
in the blue side and 32500 in the red side. For the red 
side, we used a milky flat taken using a hot blue star; for 
the blue side, we used a milky flat frame taken with the 
quartz lamp. We also observed ThAr comparison lamps 
for wavelength correction. We reduced the data with the 
MIKE pipeline written by Dan Kelson and distributed 
as part of the Carnegie Python Distribution. 

We then determine the radial velocities with the 
rvsao.xcsao task in IRAF, which cross-correlates the ob- 
ject spectrum with a set of template spectra. We use 
synthetic spectra from Munari et al. (2005) as our in- 
put velocity templates. The templates have tempera- 
tures from 3000 to 7000K, matching well with the tem- 
perature of the stars we observe. Each order was fltted 
independently, yielding a velocity and error estimate. We 
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then obtain the weighted mean of the velocity estimates. 
Typical errors for the velocity estimates are 0.5 kms^^ for 
stars that are not rotating rapidly. 

3. RESULTS 

In this section, we discuss how we determine member- 
ship and rule out background/foreground contamination 
as well as the results of optical spectroscopy and pho- 
tometry of intermediate mass stars. Table [1] lists the 
RA, Dec, photometric BVRIc magnitudes and spectral 
types of targets that we classify as members of L1641 
(with some uncertain members in the F & G range, see 
discussion in 13. 3|) . Table [2] lists objects that satisfy our 
selection criteria in 12. H and observed with OSMOS tar- 
gets, but ruled out as nonmembers. 

3.1. Spectroscopy 

We use SPTCLASS to spectral-type our targets 
(|Hernandez et al.l I2004D . a semi-automatic spectral- 
typing program. It uses empirical relations of spectral 
type and equivalent widths to classify stars. It has three 
schemes optimized for different mass ranges (K5 or later, 
late F to early K and F5 or earlier), which use different 
sets of lines. The user has to manually choose the best 
scheme for each star based on the prominent features in 
the spectrum and the consistency of several indicators. 
While SPTCLASS is insensitive to reddening and S/N 
of the spectra (as long as we have enough S/N to esti- 
mate the spectral indices), it does not take into account 
the effect of the hot continuum emission produced by 
the accretion shocks. This continuum emission makes 
the photospheric absorption lines appear weaker. SPT- 
CLASS generally assigns an earlier spectral type to veiled 
stars than their and therefore the SPTCLASS outputs 
should be considered as the earliest spectral type limits. 
While SPTCLASS utihzes spectral lines from 4000A to 
9000A, our OSMOS observations only cover a small spec- 
tral range from BIOOA to 5950A. As a result, SPTCLASS 
was unable to classify some of the K type stars, which 
we classified by eye and assigned a larger uncertainty. 
Spectra of B & A stars are also rectified and classified by 
eye as a double-check, and the results are always similar 
to the SPTCLASS output within the errors. Figure H 
shows some sample spectra. 

The earliest star found in the L1641 is a BO. 5 star right 
at the very northern edge of the field. There are three 
other Bl to B3 stars within -6.1° <Dec < -6.0°. The 
earliest star south of this range is a B4 st ar at (RA,Dec ) 
= (5:42:21.3, -08:08:00), first identified bv lRacind (fl96l . 
Among the 57 stars identified as members, there are 4 
B0-B3 stars, 14 B4 - B9 stars, 9 A stars, 15 F stars, 5 
G stars and 10 K stars. (Since we do not have radial 
velocity for our entire sample, we expect some of the A, 
F & G stars listed here are line-of-sight contamination. 
See 13X2] ') 

3.2. Photometry 

Figure |4] shows the V-I vs. V and B-V vs. B color- 
magnitude diagram of all the members listed in Table [T] 
that have magnitude information in the bands plotted. 
For stars brighter than V'^ 11 the OSMOS observa- 
tions might be saturated and the Tycho-2 magnitudes 
are used. Note that the most massive members are 



only shown in the B vs B-V CMD because we do not 
have their I band magnitude. The top panels show the 
CMD without extinction correction. The red symbols 
are objects with IR excess and the bla ck symbols are 
the non-IR excess objects, as identified bv Me geath et all 
(2012). The bottom panels show the CMD with extinc- 
tion correction, estimated by as suming intrinsic colors 
from lKenvon fc Har tmann* (*1 995l KH95) for stars earlier 
than M4 and iLegge' tt (1992) for stars M4 and later and 
a standard extinction law with Ry of 3.1 (iCardelli et al.l 
1989). Overplotted are the ZAMS (solid) and 2, 4, and 
10 Myr isochrones from lSiess et al] (j2000l ). The earliest 
type stars are essentially on the main sequence whereas 
the later type stars are still above the main sequence. 

Stars that are identified as non-members are removed 
from this figure. The background giants generally have 
large extinctions and lie above the 2 Myr isochrone 
around V-I of 0.5 - 1 and B-V of 1 in the extinction 
corrected CMD. Foreground objects are generally later 
type stars that are well-above the isochrones. 

3.3. Membership 

We expect significant contamination by non-members 
in our sample given the distance of L1641 and its large 
area on the sky. In paper I, we used indicators of youth 
such as Li absorption and Ha emission in low-resolution 
optical spectra as well as the presence of IR-excess to 
determine membership for late type stars. In the inter- 
mediate mass range (A to G), confirming membership is 
more challenging. First of all, the disk fraction is lower 
for earlier type stars so we only find a small fraction of 
members with IR-excess. Secondly, intermediate mass 
stars have stronger Ha absorption lines which can ob- 
scure Ha emission from accretion. Finally, Li depletion 
timescales are longer for intermediate m ass stars than 
for low mass stars. As demonstrated in iBriceno et al.l 
(|1997f) , G type stars in Pleiades also show Li absorption, 
and the equivalent widths are hard to distinguish from 
that of pre-main sequence stars. As a result, we have to 
use additional kinematic information to help determine 
membership, and combine this with statistical estimates 
of contamination from the Besangon model of stellar pop- 
ulations (Robin et al. 2003)as described in §3.4. 

3.3.1. Proper Motion 

One way to help eliminate foreground stars is through 
their proper motions. At the '~ 414 pc distance of Orion, 
we do not expect to detect large angular proper motions; 
stars with large proper motions are thus most likely fore- 
ground objects. In Figure [5] we plot the pro per motions 
of ou r OSMOS targets from the UCAC4 (Zac harias et al.l 
120121 ). Most of our targets, as expected, show zero or 
very small proper motions. We fit Gaussians to the 
proper motion distributions in RA and Dec. The Gaus- 
sian fit to the proper motion distribution in RA peaks 
at 0.2 mas/yr and has a width of 2.0 mas/yr; the dis- 
tribution in Dec peaks at -2.0 mas/yr with a width of 
2.7 mas/yr. 

We therefore apply a 3cr cut of 10 mas/year, or 20 
kms^^ and stars outside the circle (after taking into ac- 
count of the proper motion errors) are considered fore- 
ground objects. Note that while 90% of the objects in 
this plot has errors of < 5 mas/yr, a handful of objects 
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Figure 3. Examples of spectra from OSMOS arranged sorted by spectral types. 



Table 1 

Intermediate-mass Members of L1641 



ID 


RA 


Dec 


B 


V 


R 


I 








Spec. Type 


IR excess*^ 


RV<^ EW(Li)'= Note 




(J2000) 


(J2000) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 






(kms-i) (A) 


1 


83.31498 


-6.10956 


10.65 


10.20 


9.93 


9.67 


9.32 


9.08 


8.99 


F4.0±2.0 


N 


17.0 0.1 


2 


83.33185 


-6.28275 


13.92 


12.72 


12.08 


11.52 


10.63 


10.09 


9.93 


K1.0±3.0 


N 




3 


83.36585 


-6.04887 


9.92 


9.77 


9.71 


9.63 


9.45 


9.31 


9.29 


A4.0±2.0 


N 




4 


83.43192 


-6.40788 




13.59 




12.23 


11.21 


10.53 


10.39 


K3.0±3.0 


N 




5 


83.46720 


-6.61558 


9.47 


9.48 


9.49 


9.47 


9.31 


9.22 


9.18 


B9.0±1.5 


N 




6 


83.56561 


-6.60124 


15.21 


14.06 


13.37 


12.46 


10.59 


9.65 


8.97 


G8.0±2.0 


Y 




7 


83.61783 


-6.17349 


13.93 


12.82 


12.15 


11.47 


10.51 


9.91 


9.73 


K2.0±3.0 


N 




8 


83.67800 


-6.63059 


14.07 


12.58 


11.65 


10.72 


9.49 


8.99 


8.78 


F7.0±2.0 


N 




9 


83.69040 


-6.36113 


15.40 


14.06 


13.18 


12.27 


11.02 


10.35 


10.16 


K4.0±2.0 


N 




10 


83.70463 


-6.00634 


8.19 


8.26 






8.55 


8.60 


8.63 


B8.3±1.4 


N 





J,H,Ks photometry is from 2MASS 
^ Criteria for IR excess are defined in Gutcrmuth et al. (2009) and Megeath et al. (2012). 
From Magellan/MIKE echelle data. 
Herbig Ae star V380 Ori. 



have larger proper motion errors and are therefore com- 
patible with zero proper motion. This will make our 
effort to eliminate foreground objects less effective. 

3.3.2. Radial Velocity 

We observed 20 out of the 45 stars we classified as po- 
tential F & G members of L1641 with the MIKE echelle 
spectrograph. The radial velocities of these 20 objects 
are listed in Table [1] and [2] and shown in Figure [H Mem- 
bers of L1641 should have LSR radial velocities compat- 
ible with the gas velocity (2 - 10 kms~^, with a gradi- 
ent along north-south direction, iBallv et al.lll987l ) . Even 
though F & G stars usually have very weak Li absorp- 



tion, we are able to measure the equivalent widths at high 
spectral resolution. In order to be considered a member, 
a star should have velocities compatible with the gas (we 
chose a range of -3 - 15 kms~^, or within 5kms^^of the 
gas velocity) and also Li absorption stron ger than that 
observed in Pleiades (jBriceno et al.llT997l ). In Figure El 
we can see that most of the F& G stars observed do 
not have the cloud velocities and can be ruled out easily. 
There are 4 members that are compatible with the cloud 
velocities and show Li absorption. There are also one or 
two stars that have the cloud velocity but do not show 
Li absorption. The typical error for radial velocity de- 
termination is about 0.5 kms^^; or the size of the dots. 
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Table 2 

Non-members of L1641 Observed with OSMOS 



ID 


RA 


Dec 


B 


V 


R 


I 


J'' 






Spec. Type 


IR excess 


RV 


EW(Li) 


Note 




(J2000) 


(J2000) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 






(kms-i) 


(A) 




1 


83.24890 


-6.11734 










5.52 


4.30 


3.93 


M8.0±2.0 


N 






g 


2 


83.26690 


-6.00261 


15.31 


13.80 


12.99 


12.22 


11.00 


10.27 


10.08 


G6.0±2.0 


N 






g 


3 


83.31976 


-6.40979 


11.04 


10.45 


10.12 


9.81 


9.42 


9.12 


9.03 


F9.0±2.0 


N 






f 


4 


83.33801 


-6.41991 


12.28 


10.72 


9.89 


9.13 


7.97 


7.19 


6.97 


K3.0±3.0 


N 






g 


5 


83.36604 


-6.15616 


13.06 


11.90 


11.22 


10.56 


9.56 


8.97 


8.85 


G2.0±2.0 


N 


77.0 


0.0 


def 


6 


83.40304 


-6.38011 


10.10 


9.46 




8.79 


8.28 


8.02 


7.96 


G0.0±1.5 


N 


23.5 


0.0 


def 


7 


83.44242 


-6.46825 


13.27 


11.81 


10.99 


10.23 


9.03 


8.39 


8.18 


K0.0±3.0 


N 






g 


8 


83.48413 


-6.31374 


14.46 


13.21 


12.51 


11.84 


10.81 


10.23 


10.06 


G3.0±2.0 


N 


-18.5 


0.0 


de 


9 


83.51110 


-6.09196 


14.77 


12.88 


11.75 


10.68 


9.07 


8.13 


7.91 


F8.2±2.0 


N 


138.0 


0.0 


de 


10 


83.54718 


-6.52211 


11.81 


10.87 


10.38 


9.93 


9.23 


8.80 


8.72 


G5.0±2.0 


N 






f 



J,H,Ks photometry is from 2MASS 
^ Criteria for IR excess are defined in Gutermuth et al. (2009) and Megeath et al. (2012). 

From Magellan/MIKE echelle data. 
^ Rejected based on radial velocity incompatible with the velocity of the cloud. 

Rejected based on lack of Lithium absorption. 
^ Rejected based on large proper motion. 

^ Rejected based on position on the color-magnitude diagram. 



The larger error bars indicate fast rotating stars where 
it is more difficult to determine the radial velocity as ac- 
curat ely. However, sin ce stellar rotation decreases with 
time (|Skumanichl 1 1 9 72} ) . fast rotation is also a sign that 
these stars are young and therefore they are considered 
members. 

To summarize the result of the membership determina- 
tion, 79 out the 136 stars in our intermediate mass sam- 
ple are ruled out as non-members from kinematics and 
their position on the CMD. Only 4 out of 20, or 20%, of 
the F & G type stars observed with MIKE turn out to 
be members. In the following discussions, we will assume 
that only 20% of the stars not observed with MIKE are 
members. However, as we do not know which stars in 
our sample are non-members, we retain them in Table 1 
and Figure [71 

3.4. Completeness, Contamination and Sample 
Definition 

From here on, we combine the sample of high-to- 
intermediate mass stars characterized above and the con- 
firmed members of L1641 cataloged in paper I for our 
analysis. The sample from paper I has a total of 864 
members, including our Hcctospec and IMACS targets as 
well as the sample from Fang et al. (2009) and reference 
therein. The spectral types range from F to M4, where 
the F & G stars are mostly from the literature and not 
from our own observations. We assessed the complete- 
ness of the low-mass sample by comparing the number of 
IR-excess stars in the Megeath et al. (2012) sample and 
our sample. If we limit ourselves to objects with Av< 2 
for 1 to O.IM0, then our spectroscopic sample is about 
90% complete. Since we required Li absorption and/or 
IR-excess for membership, we expect very few contami- 
nants in the low-mass sample. We do note that the com- 
pleteness is worse for the non-excess objects because we 
are not able to target them as efhciently. In particular, 
we are missing non-IR excess objects in the G & K range 
that are too bright for 6.5m telescopes and as well as the 
lowest-mass objects because it is more difficult to ob- 
serve Li absorption in lower S/N spectra. To summarize, 
we are missing members in the low-mass sample and the 



incompleteness depends somewhat on the spectral type. 
However, in §[4l we will explain why this incompleteness 
does not affect our result as long as we use the lower limit 
of the number of low-mass stars. 

Contamination rather than completeness is the prob- 
lem for the higher mass objects, and identification of 
members is more difficult. While we can rule out fore- 
ground dwarfs and background giants to some extent 
based on their position in the color-magnitude diagram, 
but we will still have non-members in the sample that 
cannot be ruled out simply by proper motions and/or for 
which we have no radial velocity measurements. To es- 
timate the possible conta mination, we use the Besangon 
model (|Robin et al.ll2003l) to help us determine the num- 
ber and type of contaminants we would find along the 
line of sight in a 3° field toward the direction of LI 641. 

We first consider the range of O and B stars. Since 
we are trying to test the IMF variation in this mass 
range, it is essential that we have a complete sample. 
This is achieved by our survey of the brightest 2MASS 
stars, regardless of whether they have IR-excess. We 
note that because of step (3) in our target selection (see 
§ 12. ip . highly embedded O or early B stars would not 
have been found. However, evolved O or early-B stars 
would have blown away materials nearby in a manner 
similar to the B4 star near (RA, Dec) = (5:42:21.3, - 
08:08:00) and therefore would not have high extinction. 
(Note that there is another reflection nebula, powered 
by a potentially massive star further south, outside of 
the fields of our optical studies.) On the other hand, 
if there were unevolved, highly-extincted early-B stars, 
they would heat up the dust nearby and be identified as 
bright IR-excess sources by the Spitzer survey. We ex- 
amined the brightest IR-excess sources in the field and 
confirmed that they are mostly A stars, with one ex- 
ception being an FU Ori star (V883 Ori), where the IR 
luminosity is dominated by accretion. We also examined 
the protostars that have large bolometric luminosities. 
The most luminous protostar in L1641 is located at (RA, 
Dec) = (5:40:27.4,-7:27:30) i n one of the high-extin ction 
clumps. It has L ^ 490 Lg ljKrvukova et al1l2012j ). and 
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Figure 4. V-I vs V and B-V vs B CMDs of high-to-intermediate mass stars in L1641 with photometry. The dashed lines are 2, 4, and 
10 MvrlSi ess et al.l 1)20001 ') isochrones. The top panel shows CMDs without extinction correction; the bottom panel shows the same CMDs 
after correcting for extinction, assuming intrinsic colors of Kenyon & Hartmann (1995). The red symbols represent objects with IR excess. 
The open circles represent non-members and solid circles represent members. 



according to the Siess (2000) isochrones, has a mass less 
than 7 Mq. Because of the short lifetime and relative 
small number of OB stars, we expect no contamination 
in this mass range that is at the right distance to appear 
along the ZAMS at 414 pc. The Besangon model predicts 
no OB stars in the field within 1 kpc. 

For A stars, we expect our survey to be complete to at 
least Av= 2 since we targeted stars that are brighter than 
B = 15 and B-R < 2.5 in the USNOb catalog. Unfortu- 
nately, we do not have echelle observations of the A stars, 
and are therefore unable to constrain membership rate 
from radial velocities. We are left with estimating the 
level of contamination is the Besangon model, which pre- 
dicts 3 foreground A stars above the region's isochrone 
and up to 5 background giants. The actual number of 
background giants that make it into our survey depends 
on the extinction toward them and is therefore more dif- 



ficult to constrain. In Table [T] and [21 we identified 3 A 
stars as non-members and 9 A stars as members. From 
the discussion above, we therefore conclude that out of 
the 9 stars considered members, up to 5 of them can actu- 
ally be non-members. If we adopt a moderate extinction 
of Av= 2, we expect 3 non-members. 

We expect a large number of line of sight contamina- 
tion for F & G stars. In § 13.31 we determined that only 
20% of the observed F & G stars are actual members 
based on radial velocities. Our sample also suffers from 
incompleteness in the G star range because of the target 
selection criteria in our intermediate-mass and low-mass 
samples tend to miss G stars that are below the IMq line 
in the J vs. J-H CMD. 

As discussed in Paper I, our optical photometric and 
spectroscopic sample include stars in the Orion A cloud 
south of -6°, which we now define as the L1641-all sam- 
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Figure 5. Proper motion of our targets from ttie f ourtli U.S. Naval 
Obse rvatory CCD Astrograph Catalog (UCAC4; IZacliarias et al.l 
120121 1. Stars with proper motions greater than 10 mas/year, or 
201cms~^, after the proper motion error is taken into account, are 
considered foreground objects. The red and cyan symbols repre- 
sent members confirmed through IR-cxcess and radial velocity v^ith 
MIKE echelle observations, respectively. The confirmed members 
all have small proper motions, well within the 10 mas/year cutoff 
applied. 

pie. However, within this sample, the northern end is a 
relatively denser region and has four early-B stars that 
likely denote the outer regions of the ONC. We there- 
fore further divide the L1641-all sample into subsamples 
L1641-n (Dec=-6.5° - -6° and L1641-S (Dec< -6.5°). 
Note that this definition is only for the purpose of analy- 
sis in this paper and does not coincide with the small 
clusters L1641-North (RA- 84.1°, Dec- -6.3°) and 
L1641-South (RA~ 85.7°, Dec- -8.2°). 

3.5. HR diagram 

We use an HR diagram to determine the ages and 
masses of L1641 members. We use the effective temper- 
ature scale from KH95 fo r stars earlier than M 4 and ef- 
fective temperatures from lLuhman et aD (|2003l ) for stars 
M4 and later. If the Ic band magnitude is available for 
a star, we use the following equation to determine the 
bolometric luminosity: 

log(i/L0) = 0.4 [Mboi,0 - [Ic - ^/c) - DM + BCi^] 

. (2) 

where the bolometric magnitude of the sun Mboi,© is set 
to 4.75 and Aj^ is the extinction at Iq determined by 
dereddening the the observed V — Ic color. 

If the Ic band magnitude is not available (only the 
brightest stars), we use a similar equation using the V 
band magnitude to determine the bolometric luminosity; 

log(L/i0) = 0.4 [My,oi,Q-{V-Av)-DM+BCv]. (3) 

where Ay is the extinction at V determined by deredden- 
ing the B-V color. We use bolometric corrections for the 



Figure 6. LSR velocities and Li absorption data obtained from 
Echelle observations of 20 F &: G type stars in L1641. The red sym- 
bols are stars with Li absorption above the level seen in Pleiades 
llJones et al.lll99^ . A member of L1641 should have velocities 
compatible (with in 5kms~^) with the gas velocity (2-10 kms~^, 
IBallv et al.lll98^ ) and Li absorption. The typical errors in veloc- 
ity for slow-rotating stars are about 0.5kms~^, comparable to the 
size of the points. The points with error bars overplotted are fast- 
rotators, where velocity estimates are less certain. 




-10 1 2 3 4 

B-R (USNOb) 



Figure 7. USNOb B vs. B-R CMD of candidate intermediate 
mass stars selected through 2MASS. The red triangles are mem- 
bers of L1641, the black boxes are non-members of L164 (see 13.31 
for how membership is determined). The black dots are objects 
that were not observed. Note that some of the spectral types and 
membership information are obtained from Paper I. 
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corresponding magnitude band from KH95 for stars ear- 
lier than M4 and BcsseJ ()199lD for stars M4 and later. 
We assume t hat LI 641 is at the same distance as the 
ONC (414 pc: lMenten et al.ll2007l ). corresponding to dis- 
tance modulus DM = m-M = 51ogio(d/10pc)= 8.09. To 
determine the extinction, we use the V-Ic colors when 
both bands are available. For the brightest stars where 
the data is taken from the Tycho2 catalog, the B-V color 
is used to estimate the extinction. The intrins ic colors 
are from KH95 for stars earlier than M4 and iLeggettl 
(jl992|) for stars M4 and later. We then calculate the ex- 
tinct ion assuming a stan dard extinction law with Ky of 
3.1 (jCardelh et al.l[T989l) . 

Figure [8] shows the HR diagram. The black symbols 
represents stars in L1641-S and the lighter grey sym- 
bols are stars in L1641-n (see § 13.41 for th e definition 
of the populations). The 1, 3 and lOMyr 'Sie ss et aD 
(pOOO.) isochrones (Z=0.02, no overshoot) overplotted in 
red. The thick red line corresponds to the "early main- 
sequence", defined as the time when the star settles on 
the main sequence after the CN cycle has reached its 
equilibrium. In black are the evolutionary tracks for 
masses from 0.1 to 7Mq. We can estimate the ages and 
masses from the isochrones. We first select data points 
on the Siess tracks with ages less than 100 Myrs and pre- 
main sequence or main sequence evolutionary phase. We 
then use the IDL programs TRIANGULATE and TRI- 
GRID to interpolate between the evolutionary tracks to 
obtain the ages and masses. Objects that fall outside of 
the track area are not assigned masses or ages. 

3.6. Age and Mass distribution 

The top panel of Figure [3] shows the age distribution 
of the members in L1641 estimated from the HR dia- 
gram. Only objects with masses below IMq are con- 
sidered in the age estimates because the birthline effect 
tend to make th e highmass objects appear older (see e.g. 
IHartmaniill2003f ) . Both the median and the mean of the 
age distribution is around 6.5, with a standard deviation 
of 0.3 dex. Therefore, the majority of the low-mass mem- 
bers in L1641 are around 10^'^ years, consistent with the 
ages of the ONC. This is the same as the result found 
in Paper I where we used a more direct age estimate 
from the V vs. V-I CMD. This ensures that LI 641 is 
still young enough that any early B star, if ever present, 
would still be on the main sequence and seen in our sur- 
vey. There is also no evidence of SN bubbles from O 
stars. 

The bottom panel of Figure |9] shows the mass distri- 
bution of L1641 members in log-log space. Each of the 
F & G unconfirmed members is considered 0.2 star in 
the histogram. Note that the mass histogram presented 
here is not extinction limited and we did not correct for 
completeness. The dip around 1 Mq corresponds to the 
gap between the intermediate-mass sample and the low- 
mass sample. This figure demonstrates that because of 
incompleteness and extinction issues, our sample does 
not allow us to obtain the upper-mass IMF by simply 
fitting a power-low. 

4. COMPARING L1641 TO THE ONC 

We note that the selection criteria for our spectroscopic 
observations and when a star in our spectroscopic sur- 
vey is considered a member is between Paper I and this 



work and the criteria are spectral type dependent. This 
is mainly because the youth indicators themselves are 
spectral-type dependent. As a result, the completeness 
and contamination rate is also spectral- type dependent. 
As a result, our sample is not a clean extinction-limited 
complete sample and we cannot test the form of the IMF 
via the distribution of spectral types and K-band lumi- 
nosity. We can only compare the ratio of high-to-low 
mass objects. 

4.1. Comparing LI64I and ONC spectral type 
distributions 

The top panel of Figure [TU] shows the spectral type 
distribution of the extinction limited sample (Av< 2) in 
L1641 and in the ONC. The solid line shows the spec- 
tral type distribution of L1641 members with Av< 2. 
We first remove the objects that are identified as non- 
members, and then we scale the F & G stars whose mem- 
bership status is uncertain by 20%, based on the MIKE 
observations indicating that only 20% of them are actual 
members. The dashed line shows the spectral type dis- 
tribution of ONC members from H97 with Av< 2 and 
membership probability > 70%, scaled to the same num- 
ber of stars as L1641. 

The bottom panel is the same, but shows only stars in 
L1641-S. We can see that L1641, especially the southern 
part is deficient in stars earlier than B4 compared to the 
ONC. We also note that L1641 has a smaller number of 
late G and early K stars compared to the ONC. This is 
because we are missing some late G and early K stars 
from our target selection (see § 12. ip . Other than the 
missing members in this range, the low-mass end of the 
distribution look fairly similar. 

In Paper I, we found that the IMF of L1641-S is in- 
consiste nt with the stan dard IMF m o dels. Compared 
to both IChabrieil (|2005h and IKroupal (pOOlh IMFs, the 
L1641 is deficient in O and early B stars to a 3-Aa signif- 
icance level. Here we discuss whether the ratio of high 
mass to low mass stars is different in L1641 and the 
ONC. However, instead of comparing the derived masses, 
a model dependent quantity, we choose to compare the 
directly-observable spectral types. 

We then use the Fisher's exact test (jFisheii [19251 ) to 
find the significance level of the ratios between high-mass 
and low-mass stars in different regions of L1641 and the 
ONC. Since we are interested in whether the ONC has 
a higher frequency of high-mass stars, we use the one- 
sided test and do not compute the P value when the 
ONC has a lower frequency in the high-mass bin. The 
results are summarized in Tabled In general, when the 
entire L1641 population or the northern region is con- 
sidered, the frequency of high-mass stars is consistent to 
that in the ONC. In the L1641-S region, there is some 
evidence that the high-mass IMF could be different, es- 
pecially when we consider the ratio of (O to B3 stars) 
to (B4 to M4 stars), where the significance level is 0.092. 
The difference in the high-mass IMF is suggestive but not 
conclusive, mainly due to the small number of early- type 
stars in both samples. 

4.2. Comparing the LI64I and ONC K-band luminosity 

functions 

Next we compare the K-band luminosity function 
(KLF) of L1641 to that of the central regions of the ONC 
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Figure 8. HR diagr am of members of L1641. The black symbols represents st ars in L1641 - S and the lighter grey symbols are stars 
in L1641-n (see S 13.41 for the definition of the populations). The 1, 3 and 10 Mvr IS less et al.l 1)20001 ') isochrones (Z=0.02, no overshoot) 
overplotted in red. The thick red line corresponds to the "early main-sequence" , defined as the time when the star settles on the main 
sequence after the CN cycle has reached its equilibrium. In black are the evolutionary tracks for masses from 0.1 to TMq. 



Table 3 

Number of Stars in Spectral Type Bins 



Region 



O stars B0-B3 B4-B9 A stars F0-M4 



ONC (H97)=' 


2 


4 


7 


7 


437 


L164-alll'' 





4 


6 


8 


451 


L1641-n'' 





4 


3 


2 


238 


L1641-S'' 








3 


6 


213 



^ Only stars with Av< 2 and membership probability > 70% 
are considered. 

^ Only stars with Av< 2 are considered. 



(the Trapezium cluster) (jMuench et al.ll2002[ ). 

Table [5] lists the number of observed stars in the 
Trapezium and confirmed members of L1641 in each 
magnitude bin, including objects that have extinction 
higher than Ay = 2 from our optical survey as well 
as objects that are too extincted in the o ptical and are 
identified by IR-excess only. The study of 'Muenc h et al.l 
([2"002) covered the central 5'x 5'of the Trapezium and 
it is safe to assume that in the magnitude range of our 
interest (Kg brighter thanl2), the line-of-sight contami- 
nation is very small and can be ignored. In L1641, since 



the field is much larger and the density of stars much 
lower, the contaminants outnumbers the members and 
making an off-field correction can lead to large errors. 
We therefore obtain the number of all known members by 
adding the spectrally-confirmed members to the highly- 
extincted IR-excess members. The number of stars in 
each magnitude bin has varying completeness and is bi- 
ased against highly-extincted low-mass objects that do 
not show IR-excess. 

Figure [Til shows the K-band luminosity function of ex- 
tinction limited (Av< 2) L1641 sample. The top panel 



The Initial Mass Function of LI 641 



11 



Table 4 

Significance Level of Fisher's Exact Test on the Ratio of High-mass to 
Low-mass Stars in the ONC and L1641 



Compared Regions I O stars/B0-M4 0-B3/B4-M4 O-B9/A0-M4 



ONC & L1641-all'' 


0.243 


0.360 


0.314 


ONC' k. L1641n'^ 


0.421 




0.599 


ONC* & L1641s'' 


0.452 


0.092 


0.176 



Only stars with Av< 2 and membership probability > 70% are con- 
sidered. Data taken from H97. 
^ Only stars with Av< 2 are considered. 



Table 5 

Number of Stars in Ks-magnitude Bins 



Region 


Ks<7 


7 < is'^ < 8 


%<Ks<9 


9 < < 12 


ONC (Mucnch et al. (2002)) 


8 


10 


29 


397 


L1641-all 


4 


10 


23 


518 


L1641-n 


3 


2 


6 


170 


L1641-S 


1 


8 


17 


348 



Table 6 

Significance Level of Fisher's Exact Test on the Ratio of High-mass to Low-mass Stars in the 

ONC and L1641 

Compared Regions I K, <7 / 7 < Ks <12 <%/ & < Ks < 12 Ks <9 / 9 < Ks < 12 



ONC & L1641-all 


0.103 


0.118 


0.018 


ONC k Ll641-n 


0.600 


0.302 


0.050 


ONC L L1641-S 


0.034 


0.131 


0.044 



Data taken from lMuench et al.l l(2003). 

shows objects in L1641-all and the bottom panel shows 
only objects in L1641-S. 

We use the one-sided Fisher's exact test to compare the 
ratios of stars above and below the threshold magnitude. 
The hypothesis we want to test is that the L1641 has a 
lower proportion of stars in the high-mass bin compared 
to the Trapezium. We use cutoff magnitudes of Ks=7 
and 8 to ensure that we are not missing members in the 
high-mass (bright) bin. Table [5] lists the significance lev- 
els of the one-sided Fisher's test. The Trapezium popula- 
tion has a higher fraction of high-mass stars compared to 
the L1641 population as a whole even though the differ- 
ence is not very significant. The L1641-S population has 
a significantly lower fraction of high-mass stars compared 
to the Trapezium, with a P value of 3%. 

The KLF has the advantage of being minimally af- 
fected by extinction (the extinction in K band is about 
one tenth of the extinction in the V band). The KLF 
can also be affected by disk-excess. ^Mu ench et"all[200l 
constructed the K-band IR excess distribution function 
for the Trapezium, where the K band excess peaks near 
0.2 with a mean of 0.4 mag. Even though disk-excess 
can change the shape of the KLF, this excess does not 
significantly affect our number counts in Table [5] because 
of our coarse binning as well as the fact that most of the 
brightest objects are diskless. In addition, since L1641 
and the ONC has similar ages, it is reasonable to assume 
the disk excess would have similar effects on the L1641 
sample and the Trapezium sample. 



As mentioned in § 13. 4| since we can identify the higher- 
mass stars through large extinctions and we have a lower 
limit to the number of the low-mass stars, we are overes- 
timating the ratio of high-to-low mass stars and therefore 
underestimating the significance of the result. In addi- 
tion, up to half of the A stars can actually be background 
objects. If we account for these A stars, the high-mass 
bin in L1641 would have a smaller number of stars in 
the K < % and K < 9 cases. Therefore we expect the 
P value to be even smaller after accounting for A star 
contamination. 



5. DISCUSSION 
"max relation in Orion A 



5.1. Med 

iWeidner fc Kroupal (|2006l ) and iWeidner et all (|20lol ) 

compiled a list of Galactic clusters and propose a relation 
between the most massiv e star {va-ma^) in a cluster and 
the cluster mass 

(Med)- iBonnell et all ([2004f ) also pro- 
posed a similar relation based on simulations of compet- 
itive accretion simulations. In their simulations, the most 
massive star tends to form in the center of a cluster and 
gains the majority of its mass from the infalling gas onto 
the cluster, which is accompanied by newly formed low 
mass stars. Therefore, the formation of the high-mass 
stars correlates w i th hi gh stellar surface density. Sim- 
ilarly, lElmegreenI (|2004D also suggested that high mass 
stars can be formed from gravitational focused gas ac- 
cretion in high-density clouds and therefore explain the 
steep IMF in the field. The Weidner fc Kroupa (,2006l ) 
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Figure 9. Top panel: Age di stribution of obje cts with masses 
below lM0as inferred from the lSiess et al.l lj200CI 1 isochrones. We 
only consider objects with masses below lM0in the age estimates 
because the birthline effect tend to make the highmass objects 
appear older and the low-mass objects give a more reliable age 
estimate of the population. Bottom panel: Mass distribution of all 
the members. 

and iBonnell et al.l ()2004D Med - ^max relations are al- 
most identical for clusters up to a few thousand Mq, 
which suggests that, even though the authors did not 
explicitly quantify such a relation, the upper-mass IMF 
should also depend on the density. Therefore, qualita- 
tively speaking, the Meci - ^max relation is compatible 
with the environmental density dependence we find in 
Orion A. 

Quantitatively, whether L1641 follows this Meci - "camax 
relation is less clear and depends on what we define as 
a cluster. If we consider L1641-S as a whole, its most 
massive star has a smaller mass than what is expected 
from the Med - ^max relation. The L1641-S has a 
total stellar mass of '^IOOOMq, but the most massive 
star is only ~7Mq. However, if we consider the denser 



Figure 10. Top: Histogram of spectral type distribution of the 
extinction limited sample (Av< 2) in L1641-all (solid line) and 
in the ONC (dashed line). Bottom: Same as top figure, but only 
considering stars in L1641-S. The ONC population is scaled to 
match the size of the L1641 population, and the bins with only 
one star gives a negative value after scaling by a factor less than 1. 

grouping of tens of stars in L1641 individually, it would 
be consistent with the M gd - mmax relation found by 
iWeidner fc Kroupal (|2006l) . 

5.2. Challenges in studying the density dependence of 

the IMF 

Despite the apparent deficiency of high-mass stars in 
low-density regions, there are few conclusive results that 
demonstrate the IMF in low-density regions is different 
from that of clusters. Comparing the upper-mass IMF 
in a low-density star-forming region is difficult for many 
reasons, including the small number of high-mass stars 
intrinsic to the shape of the IMF, the significant line-of- 
sight contamination and the potentially different com- 
pleteness in regions of different density. 

In this work, we compared the largest nearby low- 
density and clustered star-forming regions and found 
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Figure 11. K-band luminosity function (KLF) of extinction lim- 
ited L1641 sample. Top: objects in L1641-all. Bottom: only ob- 
jects in L1641-S. In both figures, the bl ack long dashed line s rep- 
resent the KLF of the Trapezium from l IMuench et al.ll200a) . The 
solid purple lines represent the KLF of objects in L1641. 

moderately significant results. On one hand, this con- 
firms that the IMF is not vastly different in regions of dif- 
ferent density; on the other hand, the significance could 
be increased by deeper studies of the extincted popula- 
tion. We have chosen NOT to use an estimate of the 
significance that would be achieved if we made any cor- 
rection for incompleteness in the extincted low-mass pop- 
ulation to be conservative. Further study on the fainter 
members could make our result even stronger. In addi- 
tion, a complete radial velocity survey of the AFG stars 
would be useful to test IMF differences down to lower 
masses. 

6. CONCLUSIONS 

We conducted a survey of the intermediate-mass stars 
in L1641, the lower-density star- forming region south of 
the ONC, aimed at testing whether the apparent defi- 



ciency of high-mass stars in low-density regions is statis- 
tically significant. This study complements the low- mass 
survey presented in Paper I by adding 57 stars in the 
range of B4 to K4. 

In Paper I, we found that the lack of O and early B 
star s in L164 1 is in consistent with the iKroupal ()2001l ) 
and iChabriei) (|2005D IMFs. In this work, our sample 
of intermediate-mass stars improves our ability to com- 
pare directly the ratio of low-mass stars to high-mass 
stars with the ONC. In particular, we use Fisher's exact 
test to compare the spec tral type dis t ributi on of L1641 
to that of the ONC from iHillenljrandl (|1997[ ) and the K- 
band luminosity function of LI 641 directly to those in the 
centr al region of the ONC, the Trapezium (|Muench et al.l 
l2002f ). The tests indicate a probability of only 3% that 
the ONC and the southern region of L1641 (Dec<-6.5) 
were drawn from the same population, supporting the 
hypothesis that the upper mass end of the IMF is depen- 
dent on environmental density. 
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